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We consider a massless scalar field in 1+1 dimensions inside a cavity composed by 
a fixed plate, which imposes on the field a Robin BC, and an oscillating one, which 
imposes on the field a Dirichlet BC. Assuming that the plate moves for a finite 
time interval, and considering parametric resonance, we compute the total number 
of created particles inside the cavity. We generalize our results to the case of two 
parallel plates in 3+1 dimensions. 
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I. INTRODUCTION 



The dynamical Casimir effect (DCE) consists of two related phenomena: real parti- 
cle creation due to moving boundaries and radiation reaction forces on moving bound- 
aries. This effect already manifests itself for a unique moving plate and, for a non- 
relativistic motion, the frequencies of the created particles (photons in the case of the 
quantized electromagnetic field) are smaller or equal than the mechanical frequency 
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of the moving plate. Since Moore's pioneering paper [1], the DCE has been studied in 
many different situations by many authors (for a review on this subject see [2] and 
the special issue [3]). Particularly, many distinct boundary conditions (BC) have been 
considered, from the idealized Dirichlet and Neumann ones to more realistic ones. How- 
ever, the so called Robin boundary conditions (RBC), which interpolate continuously 
Dirichlet and Neumann ones, have rarely been used explicitly in the context of the 
DCE (though they have been considered by many authors in the context of the static 
Casimir effect, see for instance [4]). As far as we know, RBC appeared in the DCE 
only for the situation of one moving plate in 1+1 dimensions [5, 6]. Our purpose here 
is to consider RBC in one-dimensional cavities with one oscillating wall and in three- 
dimensional (open) cavities formed by two parallel plates with one of them oscillating 
in time. For a scalar field in 3+1 dimensions, RBC are defined by 4>\bound. = P^\bound. 
where /3 is a constant parameter with dimension of length. They interpolate continu- 
ously Dirichlet ((3 — > 0) and Neumann (j3 — > oo) BC. They appear in different areas 
of physics: from Mechanics, Electromagnetism and Quantum Mechanics to Quantum 
Field Theory, among others. These BC were used as a phenomenological model for 
penetrable surfaces [7]. In fact, for uj <C oj p , parameter j3 plays the role of the plasma 
wavelength. In Classical Mechanics, RBC may appear in a vibrating string coupled to a 
harmonic oscillator at one of its edges [5, 8]. In the context of the static Casimir effect, 
RBC lead to eigenfrequencies for the cavity modes that are roots of a transcendental 
equation. 

In the context of the Dynamical Casimir effect, Mintz et al [5] considered a massless 
scalar field in 1 + 1 dimensions under the influence of one moving boundary in a pre- 
scribed and non-relativistic motion with small amplitudes, namely, |<5g(t)| << c and 
|<5g(t)| << c/uq, where Sq(t) is the position of the moving boundary at instant t and uq 
is the dominant mechanical frequency. Using the perturbative approach of Ford and 
Vilenkin [9], the solution of the wave equation, d 2 (fi(t, x) = 0, submitted to a RBC, leads 
to a susceptibility with both real and imaginary parts, so that, SFiuj) = x(w)6Q(u)), 
with x( w ) — T^ e x{ UJ ) +i1mx(oo) ■ Recall that, for the same situation, the use of a Dirich- 
let (or Neumann) BC would lead to a purely imaginary susceptibility x D ( u ) = u 3 /6tt, 
(c = h= 1). For a typical oscillatory motion, given by Sq(t) = 5qo e~'*' //T cos(a;ot), with 
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UqT ^> 1, Mintz et al [5] showed that the dissipative force on the moving boundary 
can be enormously suppressed for (3u w 2. In a subsequent paper [6], these authors 
analyzed the particle creation phenomenon for the same situation and found that, for 
the above mentioned relation between j3 and ou , there is also an enormous suppression 
of particle creation. 



II. ONE-DIMENSIONAL CAVITIES WITH ROBIN BC 



For simplicity, we consider a one-dimensional cavity composed by a fixed plate at 
x = 0, which imposes on the massless scalar field a RBC, and a moving plate whose 
position at instant t is given by q(t), which imposes on the field a Dirichlet BC. Then, 
we must solve the wave equation with submitted to ((j>(x,t) — 7^-0(x, t) 



and t) 
namely, 



x=q(t) 
X] t), Tt(x'] t) 



■ 0. The scalar field satisfies the usual field commutation relations, 



[k(x;t),7r(x';t)} = 0, 



i5(x' — x), and (J)(x;t), (J)(x';t) 
Our anstaz for the field operators <p and fr = cj> in terms of the instantaneous basis 
introduced by Law[ll], in 1994, is 
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4>{x-t) = J2 

n 

n(x;t) = 



, ^ .. nx X]t)\h n {t) + $ n (t)\ 

t V 2k n(t) 1 J 

K{t) -„(x;t){a fl (t)-at(t)}, 



(1) 

(2) 



where the modes {u n (x, t)} of the instantaneous basis must satisfy Helmholtz equation, 
(j^2 + ^(O) u n(x;t) = 0, the boundary conditions (u n (x,t) - -f§^u n (x,t)) \ x=0 = 
and u n (x, t)\ x=q ( t ) = , and the orthonormality condition f^dxu n (x; t)u m (x; t) = S nm . 
With these properties, it follows that [cLn(t), aj„(t)] = 5 nm and [ajj(t), aj„(t)] = 
[an(t), a m (t)} = 0. The instantaneous basis can be explicitly obtained, with modes 
u n (x,t) given by 

k n (t)(x - q(t))], (3) 



/ ,s A n {t) 
u n (x; t) = — , sm 



where A n (t) = 2 



1 + 



7/<?(*) 



-1/2 



l+ 7 2fc2(t) 

dental equation sin[g(t)A; n (t)] + 7/c„(t) cos[g(t)A;„(t)] = 



and {k n (t)} are the roots of the following transcen- 



4 



Time evolution equations for a n (t) and a^if) can be found, 

k n (t) = -ik n (t)a n (t) + J2 -in (t) h (t) + A i« (*) % f (*) ' ( 4 ) 



where 



with G nm (t) := f^dx u n (x;t)u m (x;t) (an analogous equation holds for a^(t))). Re- 
lating and a* for different times, we write 

= a nm{t)a m {t ) + Y 0nm (*)4i(*o) > ( 7 ) 
rre m 

where the Bogoliubov coefficients must satisfy a nm (t ) = 5 nm and /3 nm (t ) = 0. The 
time evolution of these coefficients can be established, 

Unm{t) = -ik n (t)a nm (t) + Y Z jn (t)a jm (t) + ^ Pjm(t)i ( 8 ) 

3 3 

Pnm(t) = -ik n (t)/3 nm (t) +J2 E jn(t)Pjm(t) + ^ A jn (t)a* jm (t) . (9) 

3 3 

Previous equations may be simplified with the aid of definitions: 

K n {t) := [ dt' k n (t'); (10) 
•/ to 

^ (f\ — 6 (f\ i[K m (t)-K n (t)} . a /.\ — a -i[tf m (t)+tf n (t)] 

Consequently, the time evolution for coefficients a nm and (3 nm are 

S„m(*) = ^^„(t)^ m (t) + ^A jn (t)^ m (t); (11) 

3 3 

km(t) = J2 ^3n(t)P jm (t) + J2 hn{t)a* m {t) . (12) 
3 3 

Up to this point, our calculations are exact. However, from now on, we shall consider 
only oscillating motions with small amplitudes, so we write q(t) — qo[l + e£(i)], with 
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e 1 and £(t) given, for a typical motion, by 



sin(cu t) < t < tf 







t < or t > t f 



Expansions in powers of e (recall that all quantities get an implicit 
e-dependence through q(t)) lead to 



(13) 
(14) 



l'=\ j Jt o 



nnf„\z.*(.t-n, 



(.1=1 j Jt o 



where the superscritps mean the order of the derivative respect to e of the quantity 
in question and conditons dinm(t) = 5 nm and j3nm(t) = are satisfied. The number of 
particles created inside the cavity, with energy u n = k n , after the motion is finished is 
given by 



K(t,) = (0| 4 (t f ) K (tf) W = J2H e ^S(t f ) 

j l 

The first correction to N n (tf) occurs at order e 2 , 

m m 

For the motion in consideration, we have 

N n {tf) = ^ \Cnm(l)fnm(Uo,tf)\ 2 (eUJ tf) 2 . 



(15) 



(16) 



(17) 



where 



fnm(wo',t) :- 
Cnmin) = 



e i(u>o+K„ m )t I p—i(0JQ—K„m.)t 

(u + K nm )t 



(uJq — K nm )t 



\ A (0)A m Vk n (o)k m (Q) 



= k n (0) + k m (0). 

At this order, the total number of particles created inside the cavity is given by 
N = \Cnm{l)fnm{uo,tf)\ 2 (eu tf) 2 while the total energy of the created particles 
is given by £ = Y. k n\C nm {l)fnm{^,tf)\ 2 (euj t f ) 2 . The behavior of \f nm {^t f )\ 2 is 
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shown in Figure 1. For u tf 3> 1, it has a peak around u = K mn whose width 5 is 
proportional to l/(K nm tf) (a simple estimative gives 5 ~ 5.6/(/c nm t/)). Hence, in a first 
approximation, \ f n m(^o] tf)\ 2 behaves like a delta function, showing that whenever the 
oscillation frequency Uq equals the sum of two energy levels of the corresponding static 
cavity we have the best conditions for particle creation. 




1 ^o/ K nm 



FIG. 1: |/nm(wo; tf)\ 2 as a function of ujq/ ' n nm for ui^tf ^> 1. 



The fact that n mn is given by a sum of 2 terms means that particles are created in 
pairs. The set of values of n mn are called the resonances of the problem. Note that, 
for each value of the Robin parameter, 7, we have a different set of resonances. Figure 
2 shows how the resonances vary with 7. Since 7 varies from (Dirichlet BC) to 00 
(Neumann BC), it is convenient to make the plot against logioi^ I Qo) , instead of 7. 

For a given value of 7, the resonances are obtained by tracing a vertical line and look- 
ing at the intersections in Figure 2. The values obtained this way for log 10 (7/go) = —2 
(extreme left on the graph) are, approximately, the resonances for Dirichlet-Dirichlet 
BC since, for this case, 7 go- By the same token, the values obtained this way 
for log 10 (7/go) = 2 (extreme right on the graph) are, approximately, the resonances 
for Neumann-Dirichlet BC since, for this case, 7 ^> go- Adjacent resonances are 
equally spaced only for D-D and N-D cases. For these cases we have degenera- 
cies, which are broken in the Robin-Dirichlet case. For instance, for this last case, 
^13 = ki + k 3 7^ k 2 + k 2 = k 2 2, as can be seen in Figure 2 near log 10 (7/go) = 0. Note, 
also, the monotonic behavior of the curves with 7/go. 
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FIG. 2: Resonances (in units of ir/qo) as functions of log 10 (7/(70)- 



Figure 3 shows the number of created particles with energy k\ for different resonant 
values of the mechanical frequency as a function of log 10 (7/(70) (since particles are 
created in pairs, there are many ways of creating particles with energy k%, namely, 
1^0 = Kn, ujq = K12, etc.). For the resonance uq = «i m , we have N\ m = Cf m (j) (e/«i m t/) 2 . 

It is worth saying a few words about how the curves in Figure 3 are traced. For 
each value of 7, we compute numerically the set of corresponding resonances. Then, 
we compute Nu, _/V 12 , AT 13 , for that value of 7. We, then, take another value of 7 
and compute the new values of the resonances. Taking uq equal to the new values of 
resonances we compute again Nu, N12, Ni 3 , and so on. Hence, distinct points of a given 
curve, for instance N±, are computed with distinct values of uq, but with uq always 
equal to the first resonance (ku, which depends on 7). Note, also, the monotonic 
behavior of curves in Figure 3. 

Let us check some particular cases. For ujq = ^ and 7 = 0, which corresponds to 
the D-D case with uo = 2ki (parametric resonance with the lowest level of the static 
cavity), we have J\f\ — an( ^ £1 — ^ (l^f) ^ n agreement with Dodonov and 

Klimov [12]. For u — ^ and 7 — > 00, which corresponds to the N-D case with uq = 2k\ 
(parametric resonance with the lowest level of the static cavity, which is | the value for 
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FIG. 3: Number of created particles with energy fci(0) (in units of y-^-J f° r the first five resonant 
values of o>op as a function of log 10 (7/9o)- 



the D-D case), we have M\ — (i^f) an d — 2^ (l^f ) ' * n a g reemen t with Alves 
et al [13]. 



III. PARALLEL PLATES IN 3+1 DIMENSIONS WITH ROBIN BC 



Here, we shall generalize some of the previous results to 3+1 dimensions. We, then, 
consider a fixed plate at z = 0, which imposes on a massless scalar field a RBC and a 
moving plate, parallel to the first one, which imposes on the field a DBC. Let q(t) be 
the position of the moving plate at instant t. Operators <j> ( x ; t) and 7r (x; t) are given, 
in terms of instantaneous basis, by 



«> = E 



n=l 



d 2 kn 



7T (x; t) 



yj 2u n {k\\ , t) 

00 „ 



u n (x;t) 



ikii -x 



2n 



l n (k||,t) + h.c. 



2tt 



(k||,t) - h.c. 



(19) 



where (k\\,t) =kj + and w„ (x,t)=yj -^A n (i) sm\k n {t) (x — q(t))], with A n (t) 
and /c n (t) defined as in the 1+1 case. We shall consider the same motion as in the 1+1 
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case. The Bogoliubov coefficients are now defined by 

oo 

a n (k h t) = [a nm (k\\,t) a m (k h 0) + p nm (k h t) al (-k h 0)] . (20) 

m=0 

A perturbative solution, up to first order in e, leads to 

P nm (k h t) eM k w> = -eC nm (k\\) f nm (k h t) , (21) 

where we defined 

Cnm (k ) = ^(0)^(0) k n (0)k m (0) (22) 



U).„ 



UJ n (k\\) + UJ m 



fnm(k\\,t) = [ dt'£(t>^( fc llW fc ll)]*' (23) 



with uin(k\\) = w 2 (/cy, 0). The number of created particles in a given mode with 
k z = k n and with a parallel moment between ky and k|[+d 2 k|| is 

2 oo 

N n {k h t f ) d 2 k|| = e 2 -— 2 V \C nm (fc||) f nm ( fc|| , tf) \ d\ . (24) 

( 27r ) m=l 

The total number of created particles inside the cavity takes the form 

M{t f) = e ^ E / dk W k W \°™ ( k \\) f™ ( k btf)\ 2 > (25) 

n,m=l ® 

and the total energy is given by 

j 2 00 poo 

^( t /) = e2 ^E / dk W k \\"rn{k{)\C nm {k{)f nm {k h t f )\ 2 . (26) 
n,m=l " 

For the harmonic motion considered before, with out f ^ 1), we get 

I j, /, ^ x 1 2 TTUJ 2 tf ( 0J n (k\\) 0J m (fell) \ , - 

|/nm = [ ; , ; , 7,^ 1 S (k nm (u) - fc||) , (27) 



W n /Cm + W r , 



where y/k* m (u) + A; 2 + a/Zw^w) + fc^-w = 0. Using last result for |/ nm (k\\,t f )\ 2 , we 
obtain 

L 2 t °° 

U ^/) = e '^7 E (A n A m k n k m ) 2 Q(u-k n - k m ) (28) 

n,m=l 

and £ (tf) = Af (tf) ui/2. Figure 4 shows the behavior of the total number of created 
particles inside the plates in terms of the frequency u of the moving plate. We plot 
Af (tf) divided by e 2 u 3 L 2 tf in terms of ojao/ir. 
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0.03 — I 




FIG. 4: Total number of created particles for an open three-dimensional cavity formed by two parallel 
plates as a function of the frequency of the moving plate. 

Solid lines connected by dotted lines correspond to the DD case, while solid lines 
connected by dashed lines, to a RD case. The discontinuities occur at the resonant 
values {oj = k n + k m ). The main difference between DD and RD cases consists in the 
fact that the resonances for the former are equally spaced, while for the latter they are 
not equally spaced, as can be seen from Figure 4. Note the presence of small solid lines 
for the RD direct consequence of the degeneracy breaking that happens when 

we use RBC, as discussed previously. It is worth noting the similarity of the graph for 
the D-D case with that for the electromagnetic field inside two parallel and perfectly 
conducting plates discussed by Mundarain and Maia Neto[14]. 
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IV. FINAL COMMENTS 



In this work we considered RBC in one-dimensional cavities and in a three- 
dimensional open cavity formed by two parallel plates. Using the instantaneous basis 
method [11] we computed the number of created particles when the frequency of the 
oscillating plate was at resonance. As we showed, for one-dimensional cavities, there 
are more resonances for the RD case than for the DD or ND cases, due to the degener- 
acy breaking discussed in the text. For the same reason, there are more discontinuities 
in Figure 4 when a RBC is involved than for the case where both plates impose a DBC 
on the field. An important difference between the 1+1 and 3+1 cases treated here is 
that in the former the total number of created particles, Af (tf), is proportional to t"j, 
while in the latter, Af(tf) is proportional to tf, as shown in equations (17) and (28). 
This occurs because in the 1+1 case we have a closed cavity, while the system formed 
by two parallel plates correspond, in fact, to an open cavity. 

The possibility of suppression of the DCE [5, 6] was not investigated, since we 
considered here always resonant cavities. It would be interesting to study a massless 
scalar field in 3+1 dimensions submitted to a RBC at one moving plate and check 
if suppression of the DCE still occurs. We think that RBC, as well as more realistic 
BC, should be more investigated in the DCE, whose experimental verification seems 
imminent [15] (see also the recent proposal of experiment [16]). In this work we were 
concerned only with the regions inside the cavities, but an analysis involving also the 
outside regions, including a discussion of the dissipative force on the moving plate and 
the energy balance, can be made and will appear elsewhere. 

Acknowledgments: C.F. would like to thank CNPq and Faperj for a partial financial 
support. 
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